METHOD AND APPARATUS FOR CALIBRATION OF INDIRECT 
MEASUREMENT SYSTEMS 



Background of the Invention 

5 The present invention relates generally to instrument calibration and 

more particularly to a novel method and apparatus for calibration of indirect 
measurement and inspection systems. 

Inspection and measurement systems often depend on accurately 
sensing one or more physical parameters reflecting a process of interest. In 

10 electronics, for example, automated test equipment (ATE) may, in various 
situations, be required to accurately measure voltage, current, resistance, or 
phase. In radiation dosimetry, systems may be required to accurately 
measure total accumulated dose or radiation level. Similarly, systems may 
be required to measure length, time, mass, or any of a variety of physical 

15 parameters. It is generally desirable that the measured result be repeatable, 
portable, and accurate. Repeatability means that multiple consecutive 
readings of the same parameter on a single system give the same result 
within the required tolerance. Repeatability is primarily a function of system 
design and is beyond the scope of the present invention. Portability means 

20 that comparable results should be obtained if the same measurement is 
taken using different instances of the same system (i.e., two voltmeters of 
the same model should return similar or identical results when performing 
the same measurement). Accuracy means that the measured results should 
be a true reflection of the underlying process, without the introduction of 

25 bias. Since all measurement and inspection systems are subject, in varying 
degree, to manufacturing variations and drift, periodic calibrations of the 
measurement system are typically required to assure portability and 
accuracy. 

Frequently, calibration involves adjusting the system to obtain a 
30 correct reading on one or more known calibration samples. It is desirable 
and often required, for the calibration samples to be traceable to a primary 
standard maintained by a widely recognized agency. In the United States, 
for example, the National Institute of Standards and Technology (NIST) 
typically maintains primary calibration standards. Depending on the 
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economics of a particular situation, calibration samples may be shipped with 
the system or may be incorporated into the system itself. More commonly, 
however, small systems are returned to centralized metrology labs for 
periodic calibration, while for larger systems it may be necessary to transport 

5 the standards to the system. Such approaches are widely used and often 
suffice when the physical property of interest is directly sensed by the 
system in question, or when the relationship between the quantity sensed 
and the physical property of interest is well known. (An accurate voltage 
sensor can be combined with a known resistance to provide a current 

10 sensor, for example). 

In many cases, however, inspection and measurement systems must 
rely on measurements that only indirectly reflect the physical parameters of 
interest. The Agilent 5DX Automated X-ray Inspection (AXI) System 
manufactured by Agilent Technologies Inc. of Palo Alto, CA, for example, 

15 uses penetrating radiation, specifically x-rays, to form cross-sectional images 
of solder joints in electronic assemblies, and to automatically identify 
defective and / or unreliable joints. A key physical parameter for such 
identification is solder thickness. Unfortunately, the relationship between 
solder thickness and the gray value observed in cross-sectional or 

20 transmission images is complex and difficult to model accurately. While 
thicker joints typically appear darker, all other factors being equal, the 
relationship is complicated by several factors. Due to scatter, the gray value 
observed for a particular joint can be affected by neighboring or even distant 
portions of the assembly. Additional problems arise from the use of 

25 broadband (typically bremsstrahlung) x-ray sources in conjunction with 
monochromatic detectors. When the object under inspection comprises 
multiple materials, each material has an x-ray attenuation coefficient that 
depends on energy in a characteristic way. Monochromatic detectors sense 
only the integrated intensity reaching the detector, so it is typically not 

30 possible to unambiguously recover the thickness of the intervening 

materials. Similar problems are well known in other areas, for example, in 
various forms of quantitative medical imaging including, but not limited to, 
computed tomography (CT), positron emission tomography (PET), and 
single photon positron emission computed tomography (SPECT). 
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Calibration in such indirect-sensing systems is often expensive, time- 
consuming, and error prone. A calibration coupon, for example, comprising 
differing known values of a parameter of interest may be provided with each 
indirect-sensing system to be used as a set of calibration samples. For 
5 calibration, measurements are taken of an indirect parameter under standard 
conditions, and a fitting procedure generates a mathematical model 
describing the non-linear relation between the measured indirect parameter 
value and the actual parameter of interest. The fitting procedure must be 
repeated periodically to compensate for drift of known and unknown origin. 

10 As a specific example, consider automated X-ray inspection systems 

such as the Agilent 5DX AXI System mentioned previously. Using the prior 
. art technique, each 5DX system site is provided with a thickness coupon 
containing eight copper strips of varying thickness (including a thickness of 
zero) which intersect nine solder strips of varying thickness (including a 

15 thickness of zero), for a total of seventy-two "representative" calibration 
samples. During calibration, the gray values of each intersection point are 
measured under standard conditions, and a fitting procedure, such as one 
described in U.S. Patent 6,201 ,850 (hereby incorporated by reference for all 
that it teaches) to Heumann and owned by the present assignee of interest, 

20 constructs a mathematical model describing the non-linear relation between 
measured gray values and corresponding solder thicknesses. 

A number of disadvantages, representative of indirect measuring 
systems, are apparent in the above example. A large number of calibration 
samples may be required, resulting in large expenses. Additionally, 

25 measuring many samples accurately is time consuming, reducing system 
availability for normal use. Fabrication of the calibration samples can pose 
its own challenges. In the case of a thickness coupon, described above, 
eutectic solder is soft and easily deformed. As a result, strips whose 
thickness is controlled to the desired accuracy are not readily available 

30 commercially. This necessitates measuring the thickness of each solder 
strip using a NIST traceable instrument. However, mechanical 
measurement can alter solder strip thickness. Thickness may therefore not 
be known precisely, and may vary within a coupon or from coupon to 
coupon. Thickness measurements may suffer from portability and accuracy 
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problems as a result. Finally, the need to minimize the time that the system 
is unavailable for production use limits the sophistication of the fitting 
procedure used to describe the relationship between measured gray values 
and solder thickness. 
5 Accordingly, a calibration technique for indirect measurement and 

inspection systems is needed that avoids the aforementioned problems, yet 
is economical, automated, fast, highly repeatable, portable, and accurate. 
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Summary Of The Invention 
The present invention is a novel calibration technique for indirect 
measurement and inspection systems that solves the aforementioned 
problems of the prior art. 
5 In accordance with the invention, a reference indirect measurement 

system is calibrated based on a plurality of reference calibration samples 
having differing known values of the physical parameter of interest, for 
example using a superset calibration test coupon. Reference values for the 
reference calibration samples are obtained either by taking measurements 

10 on the reference indirect measurement system, or by simulation. A 
reference map function fitting procedure constructs a reference map 
function, f(Mc), describing the (typically non-linear) relationship between 
reference values and corresponding known parameter values. The 
reference map function fitting procedure need only be performed a little as 

15 one time on a single reference system. The resulting reference map function, 
f(Mc), is provided to each instance of one or more non-reference indirect 
measurement systems. 

Each one or more non-reference indirect measurement systems are 
also provided with a correction function fitting procedure that generates a 

20 correction function, g(M), that corrects any measurement, M, to its corrected 
value Mc to minimize differences between the local system and the 
reference system. In particular, the correction function fitting procedure 
obtains measurement values from a small number of calibration samples 
measured on the particular non-reference indirect measurement system. 

25 The correction function fitting procedure fits these obtained measurement 
values to corresponding known reference values that are associated with 
each of the calibration samples and are measured on or simulated for the 
reference indirect measurement system to generate a correction function. 
Preferably, the correction function fitting procedure fits the calibration sample 

30 points to a line/linear function, a low order polynomial, or other parametric 
function characterized by a small number of parameters to generate the 
correction function, g(M). Non-reference indirect measurement systems may 
then be corrected quickly and accurately to match the results that would 
have been obtained had the reference indirect measurement system made 
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the same measurements, thereby generating a high degree of portability and 
accuracy. 

In summary, reference map function, f(Mc), is generated once (or a 
small number of times) for a single reference system and is preferably used 

5 by all other such systems. Generation of correction function, g(M), on the 
other hand, is performed periodically or as needed on each local system to 
correct any measurement differences due to system drift, change in 
environment, etc. 

During normal operation, a given indirect measurement system 

10 recovers a raw measurement value, /W, applies the local correction function, 
g(M), to generate a corrected measurement value, Mc, and applies the 
global reference map function, f(M c )> to generate an estimate of the physical 
parameter of interest for that object of interest. The object / parameter of 
interest may then be classified into one of a plurality of classes based on the 

15 estimated value of the physical parameter for the object of interest. 
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Brief Description Of The Drawings 
A more complete appreciation of this invention, and many of the 
attendant advantages thereof, will be readily apparent as the same becomes 
better understood by reference to the following detailed description when 
5 considered in conjunction with the accompanying drawings in which like 
reference symbols indicate the same or similar components, wherein: 

FIG. 1 is a schematic block diagram of a first embodiment of a prior 
art indirect measurement system; 

FIG. 2 is a schematic block diagram of a second embodiment of a 
10 prior art indirect measurement system; 

FIG. 3 is a diagram illustrating prior art generation of a reference map 
function from a set of pairs of measured sample values and corresponding 
parameter values; 

FIG. 4 is a schematic block diagram of an indirect measurement 
15 system implemented in accordance with the invention; 

FIG. 5 is a two-dimensional graph illustrating calculation of a line 
function to arrive at a correction function, g(M), for a system being calibrated; 

FIG. 6 is a diagram illustrating generation of a reference map function 
and a correction function required for two-stage calibration in accordance 
20 with the invention; 

FIG. 7 is a block diagram of a plurality of indirect measurement 
systems illustrating the two-stage calibration process of the invention; 

FIG. 8 is a graphical plot illustrating the relationship between solder 
thickness versus gray level in an image; 
25 FIG. 9 is a graphical plot illustrating steps of solder thickness used in 

calibration samples and their corresponding expected gray level values; 

FIG. 10 is a schematic geometrical representation of a laminographic 
system to which the present invention may be applied; 

FIG. 1 1 A is a perspective view of an object having test patterns 
30 including an arrow, a circle, and a cross embedded in three different planes 
of the object; 

FIG. 11 B is a laminograph of the object of FIG. 11Awhen the focal 
point intersects the plane comprising the arrow; 
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FIG. 11 C is a laminograph of the object of FIG. 11 A when the focal 
point intersects the plane comprising the circle; 

FIG. 1 1D is a laminograph of the object of FIG. 1 1 A when the focal 
point intersects the plane comprising the cross; 
5 FIG. 1 1 E is a conventional projection radiograph of the object of FIG. 

1 1 A that produces a sharp image of each of the arrow, circle, and cross in 
the respective planes of the object; 

FIG. 12A is a schematic diagram of a typical laminographic apparatus 
in which the present invention may be employed; 
10 FIG. 12B is a top view enlargement of an inspection region shown in 

FIG. 12A; 

FIG. 12C is a perspective view of the embodiment of the 
laminographic apparatus of FIG. 12A; 

FIG. 13 is a schematic cross sectional representation of a portion of a 
15 two component assembly that may be inspected by the laminographic 
apparatus of FIG. 12A; 

FIG. 14A is a top view (in the x-y dimension) of a superset calibration 
test coupon used for calibrating a reference laminographic apparatus such 
as that of FIG. 12A to generate a reference map function f(M c )\ 
20 FIG. 14B is a side view (in the x-z dimension) of the superset 

calibration test coupon of FIG. 14A; 

FIG. 14C is a side view (in the y-z dimension) of the superset 
calibration test coupon of FIGS. 14A and 14B; 

FIG. 15 is a schematic diagram of the laminographic apparatus of 
25 FIG. 12A in which the superset test calibration coupon of FIGS. 14A-14C is 
under inspection; 

FIG. 16A is a top view (in the x-y dimension) of a two-point subset 
calibration test coupon used for performing correction function calibration of 
a local laminographic apparatus to generate a correction function g(M)\ 
30 FIG. 1 6B is a side view (in the x-z dimension) of the two-point subset 

calibration test coupon of FIG. 16A; 

FIG. 1 6C is a side view (in the y-z dimension) of the two-point subset 
calibration test coupon of FIGS. 16A and 16B; 
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FIG. 17 is a schematic diagram of a local laminographic apparatus 
which the subset test calibration coupon of FIGS. 14A-14C is under 
inspection in order to calibrate the local laminographic apparatus to the 
reference laminographic apparatus of FIG. 15; and 

FIG. 18 is a flowchart illustrating the method of measurement and 
two-stage calibration of the invention. 
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Detailed Description 
Turning now to the drawings, FIG. 1 is a simplified schematic block 
diagram of a first embodiment of a prior art indirect measurement system 10. 
As illustrated, system 10 includes one or more sensors 14 which generate a 

5 sensed signal 16 that feeds a buffer/amplifier 18. The output 20 of the 

buffer/amplifier 18 is input to a map function 22. The map function 22 may 
be implemented as program instructions stored in memory 24 and executed 
by a processor (not shown). The map function 22 may be algorithmic, 
actively calculating an estimated parameter value, P e , for each measurement 

10 value, M, or alternatively, estimated parameter values for given 

measurement values may be stored in memory, for example in a look-up 
table. The parameters and / or algorithms which need to be stored can 
reside in computer memory, in analog hardware (e.g. in settings of 
potentiometers or choice of component values), or in algorithmic form (e.g. 

15 executable code residing in firmware or software. The output 26 of the map 
function 22 is connected to a classifier circuit 28 implemented in the 
illustrative embodiment by a first comparator 34, a second comparator 36, 
and a NOR gate 38. 

In operation, a sample 12, S, possessing a physical property, P(S), of 

20 interest is observed via the one or more sensors 14. The sensed result 16 is 
conditioned/amplified by the buffer/amplifier 18 to generate one or more 
measurements 20, M. The measurement 20, M, is reflective of the physical 
property or properties, P(S), of the sample 12, S, but as described above, 
the relation may be complex and poorly understood. A stored function, f(M), 

25 approximating this relation gives an estimate 26, P e , of the physical 

parameter of interest. If desired, the estimated value 26, P e% of the physical 
parameter of interest is then passed to classifier 28 for determination of 
whether to pass or fail the sample 12 under inspection. In simple cases, for 
example, classification may be performed, as illustrated, by comparing the 

30 estimated value 26 to one or both of low and high limits, 32 Plow, and 30 
Phigh, respectively. If the map function 22, f(M), is monotonic and single- 
valued (i.e., f 1 (M) exists and is unique), the measurement 20, M, can be 
compared directly to limits M LOW = f 1 (Ptow) and Mhigh = f 1 (Phigh), for 
example. It will be understood by those skilled in the art that this description 
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should be interpreted conceptually rather than literally, and that many 
equivalent implementations are possible. 

Similarly, the measurement 20, M, may be digitized and some or all 
subsequent operations performed digitally rather than in analog circuitry, as 
5 illustrated in the system 50 in FIG. 2, which shows a second embodiment of 
a prior art indirect measurement system. Analog-to-digital (A/D) converter 
52 converts the analog measurement 20, M, to a digital signal, which is 
processed by a digital map function under the control of a processor 56. In 
this embodiment, the map function, f(M), is stored in memory 54 along with 

10 any other algorithms and data required for system operation. 

As described above, indirect measurement systems have 
conventionally been calibrated by measuring a large number, A/, of 
calibration samples, {Sj} at points, {M h P}, where i=1, ... N, which are used 
to construct a new estimate of f(M), or to improve an existing estimate, as 

15 illustrated schematically in FIG. 3. Here M\ denotes the measurement 
obtained for the t h sample, S/, having known value, P,, of the physical 
parameter of interest. 

Typically, f(M) will be chosen so that f(Mj) * P„ in addition to satisfying 
certain regularity constraints. It is often desirable for f(M) to be monotonic 

20 and easily invertible. Fitting methods, i.e., methods for estimating or 
approximating f(M), are well-known, such as described in U.S. Pat. 
6,201,850, cited in the background section, and are beyond the scope of the 
present invention. The process of: a) measuring a set of calibration 
samples, and b) constructing or revising an estimate of f(M) will be referred 

25 to hereafter as "conventional calibration". 

Within constraints dictated by economics of their marketplace, 
manufacturers typically strive to make measurement and inspection systems 
as identical as practical. Typically, they are relatively successful in this 
regard. As a result differences in the directly sensed quantities between any 

30 two systems of the same type are generally small, although not insignificant. 
They may have complex and poorly understood origins, but they will 
generally not be large. By exploiting this property, the present invention 
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constructs a two-stage calibration method for indirect measurement systems 
that avoids many of the drawbacks of the conventional approach. 

In both direct and indirect measurement systems, the difference 
between the directly measured quantities on two such systems can be 

5 typically modeled with good accuracy as a low order polynomial. This 
approach is applicable independent of the cause of the discrepancies, so 
long as the differences are small and smoothly varying. Consider, for 
example, the response of two imaging systems to a set of samples. Let the 
values measured on the first system be and that of the second system to 

10 be l 2 (j), where j denotes sample number. The error or discrepancy between 
measured values on the two systems may be denoted by 6(j) = l 2 (j) - h(j)- 
Assuming that indices j are ordered such that the UQ) are monotonic, h and 
6 can be equally well expressed as a function of for example, rather than 
of y, e.g., I2O1O)) and Because the discrepancy between systems is 

15 small, the unknown functional dependence can be expanded in power 
series, e.g.: 

5 (U) - a 0 +a?/t + a 2 li 2 + ... 
The well-known Taylor series, for example, has coefficients 
30^5(1^0) 

20 a,= dd/dh (li=0) 

a 2 = 0.5d 2 5/dh 2 (1 1=0) ... 
When the discrepancy between systems is sufficiently small and the 
functional dependence sufficiently well behaved, 5 can be accurately 
approximated using only the first two terms in the power-series expansion 

25 (i.e., by a straight line). In this case, a simple two-point calibration (i.e., 
measuring h and l 2 at two values of /, preferably corresponding to well 
separated values of l 1 suffices to correct measured values on a local system 
to match those on a reference system. Additional data points, if available 
can be used to assign confidence intervals to the fitted parameters. 

30 The present invention utilizes the calibration correction technique just 

described and applies it to indirect measurement systems to correct a each 
measurement M prior to estimating the physical parameter of interest using a 
reference map function f(Mc). 
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FIG. 4 is a schematic block diagram illustrating an improved indirect 
measurement system 100 implemented in accordance with the invention. As 
illustrated, system 100 includes one or more sensors 104 which generate a 
sensed signal 106 that feeds a buffer/amplifier 108. The output 110 of the 

5 buffer/amplifier 108 is input to a correction function 112, g(M), which may 
access a memory 114 which stores correction parameters or program code. 
The parameters and / or algorithms can reside in computer memory, in 
analog hardware (e.g. in settings of potentiometers or choice of component 
values), or in algorithmic form (e.g. executable code residing in firmware or 

10 software). The output of the correction function 1 12 is input to reference 
map function 118, f(M c ), which may access parameters or programs stored 
in memory 120. The output 122 of the reference map function 118 is, 
optionally, input to classifier 128. The classifier 128 may access a memory 
129 that stores one or more classification thresholds (e.g., a high limit 126, 

15 Phigh, and a low limit 124, Plow) for use in determining in which of a plurality 
of classes an object belongs based on the estimated parameter value of the 
object of interest. For simplicity, the details of the classifier 128 have been 
omitted. As before, this model is to be interpreted conceptually rather than 
literally, and many functionally equivalent implementations are possible. 

20 Many operations can be performed in either analog or digital fashion, for 
example. 

In operation, a sample 102, S, possessing a physical property, P(S), 
of interest is observed via the one or more sensors 104. The sensed result 
106 is conditioned/amplified by the buffer/amplifier 108 to generate one or 

25 more measurements 110, M. 

The present invention applies the correction function 112, g(M), to 
each measurement M 110 resulting in a corrected measurement 116, M c . 
Correction function 112, g(M), minimizes the small discrepancies between 
systems. Corrected measurement 116, M c , is reflective of the physical 

30 property, P(S), of the sample 102, S, but as described above, the relation 
may be complex and poorly understood. The reference map function 118, 
f(M c ), processes the corrected measurement 116, Mc, approximating this 
relation to generate an estimate 122, P e , of the physical parameter of 
interest. If desired, the estimated value 122, P e , of the physical parameter of 
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interest is used by classifier 128 to classify sample 102 into one or more 
categories (e.g. pass/fail), for example by comparison with one or more 
classification thresholds (e.g., a high limit 126, P H /gh, and a low limit 124, 
Plow)- It will be understood by those skilled in the art that this description 

5 should be interpreted conceptually rather than literally, and that many 
equivalent implementations are possible. 

The correction function 112, g(M), is preferably a simple function, 
described by a small number of parameters. Preferably, it is a linear 
function, a low-order polynomial (e.g., 4 th -order or lower), or a parametric 

10 function, such as a Gaussian or other simple functional form, characterized 
by a small number of parameters (e.g., 5 parameters or less). The reference 
map function 118, f(Mc), can, and typically will, be much more complex. 

FIG. 5 illustrates use of a simple linear function as a correction 
function, g(M) f for a system being calibrated. As illustrated, measurement 

15 values, M 1 and M 2 , are obtained for two calibration samples Si and S 2 on a 
local system of interest. Standard reference values, M 1ref and M 2 reu for the 
same two samples are obtained either by measurement on a reference 
system or by simulation, as described below. Correction function g(M) is 
then given by: 

~M 2 M Wef -M x M 2rcf 



20 M c = M 



M 2ref -M Wcf 
M 2 -M, 



M 2 -M, 

In this example, the correction function, g(M), is calculated based on only 
two measurements obtained on the local system, M 1 and M 2 , and the 
corresponding reference values, M 1ref and M 2re u from the reference system. 
The corrected value, M c , corresponding to any measurement, M, 

25 subsequently obtained during normal operation of the system of interest may 
then be found using the above equation. 

The method proposed in the present invention is not restricted to 
linear corrections. Any low order polynomial, or other parametric function 
with a small number of parameters, which can be quickly and accurately fit 

30 using a small set of measured data, and which suffices to describe the 

differences between measurements on a pair of systems, can be used in this 
manner. The present invention therefore provides a means to rapidly and 



Docket No. 10031241-1 



14 



easily correct any two or more systems, such that they return approximately 
identical values of directly measured quantities. 

In order to implement a complete mechanism for calibration, a 
reference system must be specified to which all other systems will be 

5 corrected. The reference system may be specified using various methods. 
In one embodiment, a single, potentially large, reference set of 
calibration samples (also referred to herein as a "superset calibration test 
coupon") is fabricated. This reference set of calibration samples is 
preferably a superset of and should not be confused with the small number 

10 of calibration samples (also referred to herein as a "subset calibration test 
coupon") that ship with or are incorporated into each system. Careful 
measurements of the reference calibration samples used in the superset 
calibration test coupon are taken, and the results of those measurements are 
then used to define the reference system to which all other systems are 

15 corrected. Cost of calibration is therefore significantly reduced, since only a 
single copy of the superset calibration test coupon is required, rather than 
one per customer site. Similarly, the time required to perform extremely 
careful measurements of the standards is no longer problematic, since this 
operation is performed as little as one time on a single reference system, 

20 rather than being repeated time and again on multiple systems at multiple 
sites. In principle, it is even possible to destroy the reference calibration 
samples in order to reuse the coupon materials once the measurements are 
taken. However, maintenance of the reference calibration samples may be 
useful in establishing portability between existing system models and new 

25 system models'yet to be developed. 

In another embodiment, detailed modeling and simulation replaces 
measurement to define the reference system. For an X-ray imaging system 
such as the Agilent 5DX, for example, Monte Carlo simulation using a 
coupled electron / photon transport package such as PENELOPE is 

30 appropriate. (See, F. Salvat, J.M. Fernandez-Varea, E. Acosta, and J. 
Semapu, "PENELOPE: A Code System for Monte Carlo Simulation of 
Electron and Photon transport", Issy-les-Moulineaux, France (2001) 
NEA/OECD). As in the case of careful measurement, differences between 
simulation and measurement are found to be small, so linear corrections 
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suffice. Using simulation instead of measurement avoids difficulties and 
costs associated with the fabrication of reference calibration samples. For 
example, it may allow extending the calibration to regimes in which 
fabrication of accurate calibration samples is impractical. Conversely, 

5 simulation may not be a completely faithful representation of reality. The 
resulting systematic errors can degrade accuracy. As a result, simulation is 
most appropriate for inspection and classification systems where 
repeatability and portability are more important than accuracy. 

Finally, consider construction of the reference map function, f(Mc), 

10 relating a corrected measurement, M c , to the physical parameter of interest, 
P e . The invention does not simplify this relationship, which, as described 
above, is often complex and non-linear. The invention is advantageous over 
the prior art in this regard, however, in that this construction, like 
measurement or simulation of the standard, needs to be performed only 

15 once rather than repeated periodically on each system. As a result, 

considerably greater resources (e.g., time, equipment, algorithm complexity, 
etc.) may often be invested in the fitting procedure that generates the 
reference map function, f(M c ), typically resulting in models more accurately 
reflecting the underlying relationship. 

20 Calibration of an indirect measurement system is a two-stage 

process, as shown in FIGS. 6 and 7. As illustrated in these figures, 
generation of the reference map function f(M c ) is performed as little as one 
time on a single reference system 202, while the preferably less complicated 
generation of the correction function g(M) is performed periodically on each 

25 system in the field (here illustrated by non-reference systems 21 d - 21 0 3 ), 
as required. 

FIG. 7 illustrates application of the present invention to correct a 
plurality of non-reference indirect measurement systems, 210i, 210 2 , and 
2103, based on a reference indirect measurement system, 202. While only 
30 three non-reference systems in addition to the reference system are shown 
in FIG. 7, the method is applicable to an arbitrary number of systems. As 
illustrated, reference map calibration is performed on the reference system 
202 using a superset calibration test coupon 204 that includes a number n of 
calibration samples with corresponding known values {P S u P erset jref} of the 
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physical property of interest, where j = 1 to n. A reference map function 
fitting procedure 206 fits the known values of the physical property of interest 
{PsupersetjREF} to the corresponding reference values {M S u P ersetjREF} 
measured on reference system 202, to generate or estimate a reference 
5 map function 208, f(M c ). Reference map calibration is the complex portion 
of the calibration, and the estimation of the reference map function 208, 
f(M c ), by the map function fitting procedure 206 corresponds closely to 
conventional calibration of an indirect measurement system. However, it 
differs in that estimation of the reference map function 208, f(M c ), by the 

10 map function fitting procedure 206 is performed only once on the reference 
system 202. This reference map function 208, f(M c ), is then utilized for all 
other systems, 210i, 210 2 , and 210 3 , referred to herein as "non-reference 
indirect measurement systems". 

Non-reference indirect measurement systems 21 0i 21 0 2 , and 21 0 3 

15 are shipped with (or otherwise obtain) the reference map function 208, a 
correction function fitting procedure 214, and respective subset calibration 
test coupons, 212 if 212 2 , and 212 3 . Preferably, non-reference indirect 
measurement systems 21 d 21 0 2 , and 210 3 are also shipped with reference 
values {MsubsetjREF}, i = 1 to m, expected for measurements on the subset 

20 coupon. Thus, in the linear calibration example of FIG, 5, as few as 2 
reference values, M SU bseuREF, M SU bset_2REF could be provided to the non- 
reference indirect measurement systems 21 0i 21 0 2 , and 21 0 3 . In the 
alternative, as described hereinafter, the reference values supplied may be 
computed by extrapolation or interpolation if the components on the subset 

25 coupon are not also found on the superset coupon. 

It should further be noted that the number of subset calibration 
samples may even be as few as one sample. In this case, the correction 
function fitting procedure would be as simple as finding an additive 
correction (i.e., an offset) to fit the measured value to the known reference 

30 value of the sample. 

Correction function calibration of each system 21 0i 21 0 2 , and 21 0 3 is 
respectively performed prior to first normal use and periodically thereafter as 
required to generate respective correction functions, 21 61, 21 6 2 , and 21 6 3 
specific to the corresponding system to bring the measurements obtained on 
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the respective systems, 21 Ch, 210 2 , and 210 3 into agreement with those 
which would have been obtained on reference system 202 given the same 
samples. To accomplish this, a reduced set of samples comprising the 
respective subset calibration test coupons, 212i, 212 2 , and 212 3 are 

5 measured by the respective systems, 21d, 210 2 , and 210 3 and fitted by the 
respective correction function fitting procedure 214 to generate respective 
correction functions, 216i, 216 2 , and 216 3 specific to their respective 
systems 2^0^ 21 0 2 , and 21 0 3 . 

As just described, the periodic correction function calibration therefore 

10 consists of conducting a greatly reduced number of measurements and 
fitting the measurements to construct a simple correction function, g(M). 
When it suffices to implement the correction function, g(M), using a linear 
function, for example, the reduced measurement set can consist of as few as 
a pair of data points, as illustrated in FIG. 5. By comparing measurements 

15 obtained on these samples to the corresponding measurements on the 

reference system, the correction function g(M) can be uniquely defined for a 
given system. Typically the samples utilized for periodic correction function 
calibration will be a subset of the larger group utilized for estimation of 
reference function f(M c ), although this is not strictly required. When the 

20 periodic calibration samples are not a subset of the larger set of calibration 
samples used to generate the reference map function, f(M c ), on the 
reference system, interpolation or extrapolation in f 1 (P) can be used to 
identify behavior of the reference system at the measured points. 

In summary, the fact that the measurement errors between indirect 

25 measuring systems are typically small is utilized in the present invention to 
provide an improved two-stage calibration technique. The improved 
technique reduces system cost, while improving portability and system 
availability. When required, accuracy may also be improved, since the 
complex portion of the calibration need be performed only once under 

30 laboratory conditions. 

Turning now to a specific illustrative embodiment, the invention will be 
described in the context of an Automated X-ray Inspection (AXI) System, for 
example, the Agilent 5DX AXI System, manufactured by Agilent 
Technologies, Inc., of Palo Alto, CA. In this illustrative embodiment, the 
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parameter of interest, P, is the thickness, f, of solder, and the sensed or 
measured value, /W, is an image gray level or X-ray intensity level. 

While the following description is presented in terms of a two 
component assembly or portion of an assembly comprising a layer of solder 
5 and a layer of copper, it is to be understood that the present invention also 
applies to any two component assembly or portion of an assembly. It is to be 
further understood that the present invention applies equally to assemblies 
with three or more components. For example, when the third component is 
present and unchanging (e.g., the G10 substrate of a printed circuit 

10 assembly, its effect is simply to alter the effective source intensity spectrum. 
As such, it is not explicitly treated in the following description. The invention 
may also be extended to systems characterized by three or more variable 
components as well. Furthermore, the two components need not be in 
distinct layers but may be intermixed. One skilled in the art will recognize 

15 that the terms "gray level" and "intensity", as used throughout, are closely 
related. "Gray level" refers to a numeric value assigned to represent the 
brightness (or darkness) of a pixel in an image. Gray levels can be scaled 
arbitrarily, but a range of 0 to 255 is common. Similarly "intensity" as used in 
this example refers to energy per unit area or brightness of the X-ray beam. 

20 Modern X-ray imaging systems respond predictably and often linearly to 
given X-ray intensity. A specific image gray level is therefore functionally 
related to a corresponding X-ray intensity. While the two quantities are 
related, typically in a known way, they should not be considered identical, 
however. 

25 In an X-ray image of solder material, typically a combination of lead 

and tin, a relationship exists between X-ray intensity or gray value and the 
thicknesses of the solder imaged. FIG. 8 illustrates an example of this 
general relationship. In this example, it is seen that the image of a thin 
section of solder will have a gray level that is less than the gray level of the 

30 image of a thicker section of solder, i.e. the image of the thin section will 
appear to be a lighter shade of gray than the image of the thicker section. 
This presumes an imaging system such as the 5DX which presents X-ray 
images as positive rather than negative. Additionally, in the following 
description, the gray scale ranges from zero to a maximum value G M ax, with 
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lower values correspond to lighter pixels and higher X-ray intensity, while 
values near the maximum value correspond to darker shades of gray and 
lower X-ray intensity. It is to be understood that other conventions for 
representing the gray scale may also be used. For example, lower values 
5 may be selected to correspond to the darker shades of gray (black) and the 
values near the maximum value may be selected to correspond to lighter 
shades of gray (white). 

The relationship between solder thickness, f, and image gray level, G, 
may be calibrated using multiple but differing values of solder thickness. For 

10 illustrative purposes, FIG. 9 is a graph illustrating the relationship between 
gray value and solder thickness, with thicknesses ranging from 0.001 inch to 
0.010 inch in increments of 0.001 inch. Since the solder thickness, t, at each 
of the steps is known, the corresponding gray level values, G, may be 
compared to the gray level values, G, of other X-ray images of solder 

15 material where the thickness values, t, are not known in order to determine 
the unknown thickness values, t. 

In the case of circuit board assemblies, the solder is attached to a 
circuit board. Thus, gray scales displayed in the X-ray images include 
contributions from the solder as well as from the material comprising the 

20 circuit board. Typically the circuit board substrate is a plastic or resin type 
material and may further include ground planes and circuit traces made of a 
conducting material such as copper. In these cases, determination of the 
solder thickness is complicated by the presence of the circuit board and 
associated materials which contribute to background shading in the X-ray 

25 images. Background shading correction techniques for removing or 
minimizing the effects of such shading are described below. 

FIG. 10 shows a schematic representation 300 of a typical 
laminographic geometry which may be used with the present invention. An 
object 310 under examination, for example, a circuit board, is held in a 

30 stationary position with respect to a source of X-rays 320 and an X-ray 
detector 330. Synchronous rotation of the X-ray source 320 and detector 
330 about a common axis 340 causes an X-ray image of the plane 360 
within the object 310 to be formed on the detector 330. The image plane 360 
is substantially parallel to the planes 362 and 364 defined by the rotation of 
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the source 320 and detector 330, respectively. The image plane 360 is 
located at the intersection 370 of a central ray 350 from the X-ray source 320 
and the common axis of rotation 340. This point of intersection 370 acts as a 
fulcrum for the central ray 350, thus causing an in-focus cross-sectional X- 

5 ray image of the object 310 at the plane 360 to be formed on detector 330 as 
the source and detector synchronously rotate about the intersection point 
370. Structure within the object 310 which lies outside of plane 360 forms a 
blurred X-ray image on detector 330. 

In the laminographic geometry shown in FIG. 10, the axis of rotation 

10 of the radiation source 320 and the axis of rotation of the detector 330 are 
coaxial. However, it is not necessary that these axes of rotation of the 
radiation source 320 and the detector 330 be coaxial. The conditions of 
laminography are satisfied and a cross-sectional image of the layer 360 will 
be produced as long as the planes of rotation 362 and 364 are mutually 

15 parallel, and the axes of rotation of the source 320 and the detector 330 are 
mutually parallel and fixed in relationship to each other. It is to be understood 
that the present invention is not limited to any specific laminographic 
configuration. One skilled in the art will recognize that there are numerous 
alternative configurations for generating laminographic images which may 

20 also be used. Furthermore, the present invention is not limited to cross- 
sectional images of a two component assembly, but may be practiced with 
any type of X-ray image of the assembly, including but not limited to 
laminographic images, CT images, shadowgraph images, etc. 

FIGS. 1 1 A-1 1 E show laminographs produced by the above described 

25 laminographic technique. The object 310 shown in FIG. 1 1 A has test 

patterns in the shape of an arrow 381 , a circle 382 and cross 383 embedded 
within the object 310 in three different planes 360a, 360b and 360c, 
respectively. 

FIG. 1 1 B shows a typical laminograph of object 310 formed on 
30 detector 330 when the point of intersection 370 lies in plane 360a of FIG. 

11 A. The image 390 of arrow 381 is in sharp focus, while the images of other 
features within the object 310, such as the circle 382 and cross 383 form a 
blurred region 391 which does not greatly obscure the arrow image 390. 
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Similarly, when the point of intersection 370 lies in plane 360b, the 
image 392 of the circle 382 is in sharp focus as seen in FIG. 1 1 C. The arrow 
381 and cross 383 form a blurred region 393. FIG. 1 1 D shows a sharp 
image 394 formed of the cross 383 when the point of intersection 370 lies in 

5 plane 360c. The arrow 381 and circle 382 form blurred region 395. 

For comparison, FIG. 1 1E shows an X-ray shadow image 396 of 
object 310 formed by conventional projection radiography techniques. This 
technique produces a sharp image of each of the arrow 381, circle 382 and 
cross 383, respectively, which overlap one another. FIG. 1 1 E vividly 

10 illustrates how multiple characteristics contained within the object 310 may 
create multiple overshadowing features in the X-ray image which obscure 
individual features of the image. 

FIG. 12A illustrates a schematic diagram of a typical laminographic 
apparatus 400 usable with the present invention. In this configuration, an 

15 object under inspection is a printed circuit board 410 having multiple 
electronic components 412 mounted on the board 410 and electrically 
interconnected via electrical connections 414 (See FIG. 12B). Typically, the 
electrical connections 414 are formed of eutectic solder. However, various 
other techniques for making the electrical connections 414 are well known in 

20 the art and even though the invention will be described in terms of solder 
joints, it will be understood that other types of electrical connections 414 
including, but not limited to, conductive epoxy, mechanical, and non-eutectic 
bonds may be inspected utilizing the invention. FIG. 12B, which is a top view 
enlargement of a region 483 of the circuit board 410, more clearly shows the 

25 components 412 and solder joints 414, which in the illustrative embodiment 
are ball grid array joints. 

The laminographic apparatus acquires cross-sectional images of the 
solder joints 414 using the previously described laminographic method or 
other methods capable of producing equivalent cross-sectional images. The 

30 cross-sectional images of the solder joints 414 are automatically evaluated 
to determine their quality and physical characteristics, including, for example, 
solder thickness, t. Based on the evaluation, a report of the solder joint 
quality and physical characteristics is presented to the user. 
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The laminographic apparatus, as shown in FIG. 12A, comprises an X- 
ray tube 402 which is positioned adjacent printed circuit assembly 410. 
Printed circuit assembly 410 is optionally supported by a fixture 420. Printed 
circuit assembly 410 or fixture 420, if present, is attached to a positioning 

5 table 430 which is capable of moving the fixture 420 and board 410 along 
three mutually perpendicular axes, X, Y and Z. A rotating X-ray detector 440 
comprising a fluorescent screen 450, a first mirror 452, a second mirror 454 
and a turntable 456 is positioned adjacent the circuit board 410 on the side 
opposite the X-ray tube 402. A camera 458 is positioned opposite mirror 452 

10 for viewing images reflected into the mirrors 452, 454 from fluorescent 
screen 450. A feedback system 460 has an input connection 462 from a 
sensor 463 which detects the angular position of the turntable 456 and an 
output connection 464 to X and Y deflection coils 481 on X-ray tube 402. A 
position encoder 465 is attached to turntable 456. The position sensor 463 is 

15 mounted adjacent encoder 465 in a fixed position relative to the axis of 
rotation 440. The camera 458 is connected to a computer 470 via an input 
line 476. The computer 470 includes the capability to perform high speed 
image analysis. An output line 478 from the computer 470 connects the 
computer to positioning table 430. 

20 The master computer 470 is connected to access a memory 496 that 

stores software routines in the form of program instructions, parameters, and 
data. In the illustrative embodiment, correction function 498 and reference 
map function 499 may be stored in memory 496 as procedures executed by 
a computer processor for computing the respective functions, as lookup 

25 tables, as parametric or non-parametric function representations, or in any 
other convenient form. 

A perspective view of the laminographic apparatus is shown in FIG. 
12C. In addition to the X-ray tube 402, circuit board 410, fluorescent screen 
450, turntable 456, camera 458, positioning table 430 and computer 470 

30 shown in FIG. 12A, a granite support table 490, a load/unload port 492 and 
an operator station 494 are shown. The granite table 490 provides a rigid, 
vibration free platform for structurally integrating the major functional 
elements of the laminographic apparatus, including but not limited to the X- 
ray tube 402, positioning table 430 and turntable 456. The load/unload port 
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492 provides a means for inserting and removing circuit boards 410 from the 
machine. The operator station 494 provides an input/output capability for 
controlling the functions of the laminographic apparatus as well as for 
communication of inspection data to an operator. 

5 In operation of the laminographic apparatus as shown in FIGS. 12A 

and 12C, high resolution, cross-sectional X-ray images of the solder joints 
414 connecting components 412 on circuit board 410 are acquired using the 
X-ray laminographic method previously described in reference to FIGS. 10 
and 11A-11E. Specifically, X-ray tube 402, as shown in FIG. 12A, comprises 

10 a rotating electron beam spot 485 which produces a rotating source 480 of 
X-rays 482. The X-ray beam 482 illuminates a region 483 of circuit board 
410 including the solder joints 414 located within region 483. X-rays 484 
which penetrate the solder joints 414, components 412 and board 410 are 
intercepted by the rotating fluorescent screen 450. 

15 Dynamic alignment of the position of the X-ray source 480 with the 

position of rotating X-ray detector 440 is precisely controlled by feedback 
system 460. The feedback system correlates the position of the rotating 
turntable 456 with calibrated X and Y deflection values stored in a look-up 
table (LUT). Drive signals proportional to the calibrated X and Y deflection 

20 values are transmitted to the steering coils 481 on the X-ray tube 402. In 

response to these drive signals, steering coils 481 deflect electron beam 485 
to locations on a target anode 487 such that the position of the X-ray source 
spot 480 rotates in synchronization with the rotation of detector 440 in the 
manner previously discussed in connection with FIG. 10. 

25 X-rays 484 which penetrate the board 410 and strike fluorescent 

screen 450 are converted to visible light 486, thus creating a visible image of 
a single plane within the region 483 of the circuit board 410. The visible light 
486 is reflected by mirrors 452 and 454 into camera 458. Camera 458 
typically comprises a low light level closed circuit TV (CCTV) camera which 

30 transmits electronic video signals corresponding to the X-ray and visible 
images to the computer 470 via line 476. 

The computer 470 sends measurement values 493, /W, to the 
correction function 498, g(M) t in memory 496 which returns corresponding 
corrected measurement values M c 494. The corrected measurement values 
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494, M c , are submitted to the reference map function 499, f(M c ) in memory, 
which returns corresponding estimated solder thickness values 495, t e . A 
classification function. 497, determines whether to pass or fail the 
corresponding solder joint(s). In the simplest form of classifier, classification 

5 consists simply of comparing estimated thickness values 495, t ei to upper 
and lower threshold limits stored in memory (not shown). 

Computer 470 includes one or more processors, one or more 
memories and various input and output devices including but not limited to 
monitors, disk drives, printers and keyboards. It is to be understood that the 

10 image analysis methods of the present invention including the correction 
function 498 and reference map function 499 may be implemented in a 
variety of ways by one skilled the art. Additionally, it is to be understood that 
the term "image" is not limited to formats which may be viewed visually, but 
may also include digital or analog representations which may be acquired, 

15 stored and analyzed by the computer. 

Computer 470 also controls the movement of positioning table 430 
and thus circuit board 410 so that different regions of circuit board 410 
containing solder joints of interest may be automatically positioned within 
inspection region 483. 

20 FIG. 13 shows a schematic cross sectional representation of a portion 

of a two component assembly 500 comprising a first material 510 (e.g., 
solder) in combination with a second material 520 (e.g., copper, plastic, etc). 
As shown in FIG. 13, X-rays 530 having a incident intensity l 0t are directed 
upon the assembly 500 from a first side and encounter regions of the 

25 assembly 500 which include the first material 510 having a thickness ti in 
combination with the second material 520 having a thickness t 2l and other 
regions of the assembly 500 which include only the second material 520. In 
regions where the X-rays have passed through only the second material 
520, the incident intensity l 0 is attenuated to an intensity Similarly, in 

30 regions where the X-rays have passed through both the first material 510 
and the second material 520, the incident intensity l 0 is attenuated to an 
intensity / 2 . Ignoring scattering, the absorption of monochromatic X-rays in 
the region including only the second material 520 is governed by the 
following relation: 
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where 0C2 is the linear X-ray attenuation coefficient for the second material 
520. The absorption of monochromatic X-rays in the region including both 
the first material 510 and the second material 520 is governed by the 
5 following relation: 

where cci is the linear X-ray attenuation coefficient for the first material 510. 
FIG. 13 illustrates the X-rays 530 passing through the assembly 500 in a 
direction which is perpendicular to the first and second layers 510 and 520, 

10 thus, U and t 2 represent the thicknesses of the first and second layers 510 
and 520, respectively. In the event the X-rays pass through the assembly at 
some other angle, t 1 and t 2 represent the distances the X-rays have traveled 
through the first and second layers 510 and 520, respectively. FIGS. 14A- 
14C respectively illustrate a top view (in the x-y dimension), a side view (in 

15 the x-z dimension), and a side view (in the y-z dimension) of a superset 
calibration test coupon 550 used in calibrating a reference AXI system to 
generate the reference map function f(M 0 ). The superset calibration test 
coupon comprises eight copper strips 551, 552, 553, 554, 555, 556, 557, 
558 of varying thickness (CO, C1, C2, C3, C4, 5, C6, and C7, where CO is a 

20 minimum thickness (which may even include zero thickness) and C7 is a 
maximum thickness) intersecting nine solder strips 561, 562, 563, 564, 565, 
566, 567, 568, 569 of varying thickness (SO, S1, S2, S3, S4, S5, S6, S7, and 
S8, where SO is a minimum thickness (which may even include zero 
thickness) and S8 is a maximum thickness), for a total of seventy-two 

25 intersections which comprise the "representative" calibration samples. The 
strips are assembled into a printed circuit board assembly comprised of 
materials 570 widely used in printed circuit boards (e.g., a G10 substrate). 

Various methods exist for implementing the reference map function 
fitting procedure 472, and will of course depend on the particular subject 

30 application or measurement system. In the case of the Agilent 5DX AXI 
System, the fitting procedure for generating f(M c ) is preferably performed 
according to one or more of several fitting procedure embodiments 
described in U.S. Patent 6,201,850 (incorporated herein by reference for all 



Docket No. 10031241-1 



that it teaches) to Heumann and owned by the assignee of interest. 
Generally, the fitting procedure 472 involves a non-linear shading correction 
technique based on the assumptions that 1) plots of delta gray level (y=AG, 
the difference between foreground and background gray levels) at particular 

5 solder thickness on the y-axis vs. background gray level on the x-axis may 
be approximated by points located on a left branch of a series of hyperbolic 
curves having two common parameters including a common x-axis value at 
which each hyperbola assumes its minimum value of y (i.e., AG), and a 
common x-axis intercept (x=xo), and 2) at a "nominal" or reference 

10 background gray level, the delta gray level due to solder (y 0 ) vs. solder 

thickness (t), function may be approximated by a fitted curve of known form, 
e.g., a sum of exponentials. The procedure involves fitting multiple sets of 
calibration data to hyperbolic functions of the form: 

y = AG = JO-a) 2 +6 2 + c (3) 
15 where a is the common (non-physical) x-axis value at which each 

hyperbola has a minimum value, and all of the hyperbolae share a common 

x-axis intercept at x 0 = a-4c 2 -b 2 . Here, y represents the delta-gray level 

values, x represents the background gray level, and each hyberbola 
corresponds to a particular solder thickness. Once a set of calibration 

20 curves are fitted to the sample data, these curves are used as the reference 
map f(M c ) to map measured gray values to solder thickness during normal 
operation of the AXI system. 

To generate the reference map function, f(Mc), the superset 
calibration test coupon printed circuit board 550 is inserted into a reference 

25 AXI system 400a, as illustrated in FIG. 15, which measures the gray values 
of X-ray images of each intersection (marked with an "x" on the test coupon 
550) under standard conditions. The reference AXI system 400a is a system 
such as the system 400 of FIG. 12A that is designated as the reference to 
which all other identical-by-design systems will be calibrated. Accordingly, in 

30 FIG. 15, all identical components are shown with the same reference signs 
as those shown in FIG. 12A. A reference map function fitting procedure 472, 
executing on a remote computer system 401 , generates a mathematical 
model, i.e., the reference map function f(Mc), describing the non-linear 
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relation between gray value measured on the reference system, M, and 
known solder thickness, t. 

In the alternative, as discussed previously, detailed modeling and 
simulation may replace obtaining superset calibration test coupon sample 

5 measurements to define the reference system, for example, using the 
PENELOPE Monte Carlo simulation package described previously. 

Once a reference map function 499, f(M c )> is generated, it may be 
stored in the AXI system memory 496 as either a software routine 
comprising programmable instructions for execution by the AXI system 

10 hardware, as a lookup table 473 comprising pairs of gray values and 

corresponding solder thicknesses, {M h t}, where i=1 to N, as a parametric or 
non-parametric function representation, or in any convenient form. When 
the reference map function 499, f(M c ), is complex, as it typically will be in 
AXI systems, it will often be preferable to generate a lookup table 473 which 

15 may be accessed and interpolated more efficiently than f(M c ) can be 
computed. 

As discussed previously, once initially calibrated, an AXI system 
should be periodically recalibrated (i.e., correction function calibration should 
be performed) to ensure portability and accuracy. However, as also 
20 previously discussed, recalibration using a reference superset test coupon 
550 and the fitting procedure 472 just described is a time-consuming and 
costly process. Accordingly, the present invention introduces a simplified 
correction calibration procedure for calibration any time after the initial 
factory calibration. 

25 In the illustrative embodiment, the correction function fitting procedure 

involves the use of a simplified test coupon, for example a two-point test 
coupon 580 illustrated in FIGS. 16A, 16B, and 16C, and a linear calculation 
function 498 for generating correction function g(M) that calculates a straight 
line function using measurement data from each of the two samples or 

30 intersection points (marked with an "x") on the coupon. 

As shown in FIGS. 16A, 16B, and 16C, the two-point test coupon 580 
comprises two copper rectangles, dots, or other shapes, 581 a, 581 b, of 
predefined thickness (e.g., C7) overlaid by two corresponding solder 
rectangles, dots, or other shapes 582a, 582b, of varying thicknesses (e.g., 
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SO and S9), for a total of two "representative" calibration samples (located at 
the two intersections marked with an "x" in FIG. 16A). Preferably, the 
thickness of the two solder samples 582a, 582b are widely separated, 
spanning the desired operating range of the systems to be calibrated, where 
5 one may even have zero thickness. The copper/solder samples 581a/582a, 
581b/582b are assembled into a printed circuit board assembly 580 
comprised of materials 584 widely used in printed circuit boards (e.g., a G10 
substrate). 

The correction function fitting procedure 492 that generates the 

10 correction function 498, g(M), can simply determine the equation of the line 
mapping the first and second measured values to the corresponding 
reference values. The correction function fitting procedure 492 can be 
implemented either as a software routine residing in system memory 496 
comprising program instructions executed by the master computer 470 or as 

15 a software routine comprising program instructions residing in the memory of 
a remote computer 403, or may be implemented in hardware. Preferably, 
since only two well-separated points are required, thicknesses can be 
chosen to minimize the difficulties in production of the calibration samples. 
For example, one can choose to use air (zero thickness) as one sample, and 

20 thick copper plus thick solder as the other sample. Measured gray values at 
high thickness are relatively insensitive to minor changes in thickness. 
Additional data points, if implemented on the subset calibration test coupon, 
can be used to assign confidence intervals to the parameters of the fitted 
line. It has been found that a simple, linear correction such as that just 

25 described, accounts for greater than 99.5% of the variation between Agilent 
5DX systems, despite the fact that these are large, complex systems. Thus, 
by measuring only two calibration samples, rather than the large number 
used in the superset calibration test coupon to determine the reference map 
function f(Mc), such systems can be brought into excellent agreement with 

30 each other more quickly, more economically, and sometimes more 
accurately than with conventional calibration methods. 

As stated previously, the method proposed in the present invention is 
not restricted linear corrections. Any low order polynomial, or other 
parametric function with a small number of parameters, which can be quickly 
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and accurately fit using a small set of measured data, and which suffices to 
describe the differences between measurements on a pair of systems, can 
be used in this manner. 

To generate the correction function, g(M), the subset calibration test 

5 coupon printed circuit board 580 is inserted into a local AXI system 400b 
(which is generally a different system from, but identical in design as, the 
reference AXI system 400a), as illustrated in FIG. 17, which measures the 
gray values of X-ray images of each of the two test points, or intersections 
(marked with an V on the test coupon 580) on the two-point test coupon 

10 580 under standard conditions. The local AXI system 400b is a system such 
as the system 400 of FIG. 12A that is to be calibrated to the reference 
system 400a of FIG. 15. Its main functional components are identical by 
design to the reference system of FIG. 15; accordingly all identical 
components are shown with the same reference signs as those shown in 

15 FIG. 15. 

The AXI system software includes a correction function fitting 
procedure 492, which receives the calibration sample data 491, (M, t), for 
each subset calibration test coupon calibration sample, and generates a 
correction function 498, g(M) t which can take any gray level value of M and 

20 correct it to its calibrated value Mc. The correction function fitting procedure 
492 fits the calibration sample points and associated reference values to a 
line to generate g(M). Because g(M) is a simple linear equation, a low order 
polynomial, or parametric function with a small number of parameters, the 
correction function 498, g(M), can be and preferably is efficiently 

25 implemented in software rather than stored as a lookup table. However, 
lookup table or alternative representations for g(M) can also be used. 

In summary, the present invention provides a system and method for 
determining an estimated value of a parameter of interest of an object. FIG. 
18 is a flowchart illustrating the measurement method 600 of the invention, in 

30 which a raw measurement that is indirectly representative of a parameter of 
interest of the object is obtained (step 601 ), the raw measurement is 
corrected to a corrected measurement to minimize measurement differences 
between the indirect measurement system and a reference indirect 
measurement system (step 602), and a value of the parameter of interest of 
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the object is estimated based on the corrected measurement (step 603). 
Preferably, the physical parameter value is used to classify the object into 
one of a plurality of classes (step 604). 

In the preferred embodiment, as discussed previously, the step for 

5 correcting the raw measurement to the corrected measurement utilizes a 
correction function in the form of either a linear equation, a low-order 
polynomial , or a parametric function characterized by a small number of 
parameters (such as a Gaussian or other simple functional form) (step 607) 
to determine the corrected measurement. 

10 The correction function is determined (step 606) based on 

measurements of one or more calibration samples obtained (step 605) using 
the indirect measurement system and known reference values associated 
with each of the one or more calibration samples. 

The reference map function is preferably determined (step 609) based 

15 on measurements of a plurality of calibration samples obtained (step 608) 
using the reference indirect measurement system and corresponding known 
values of the physical property of interest associated with each of said 
plurality of calibration samples. Although this preferred embodiment of the 
present invention has been disclosed for illustrative purposes, those skilled 

20 in the art will appreciate that various modifications, additions and 

substitutions are possible, without departing from the scope and spirit of the 
invention as disclosed in the accompanying claims. The laminographic 
geometry and apparatus shown and described with reference to FIGS. 10-12 
are typical of that which may be used in conjunction with the present 

25 invention. However, specific details of these systems are not critical to the 
practice of the present invention, which addresses a two-stage method for 
calibration of indirect measurement systems in general. For example, the 
number of computers and delegation of tasks to specific computers may vary 
considerably from system to system as may the specific details of the 

30 parameter of interest, the indirect sensor mechanism, etc. The present 
invention is applicable to any type of system which derives a parameter of 
interest from an indirect measurement (i.e., by measuring a different 
parameter and mapping the measured parameter to the parameter of 
interest), and is not limited to AXI systems which have been described 
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herein for purposes of illustration. It is also possible that other benefits or 
uses of the currently disclosed invention will become apparent over time. 
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